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ABSTRACT 

A pulsar wind nebula inside a supernova remnant provides a unique insight into 
the properties of the central neutron star, the relativistic wind powered by its loss 
of rotational energy, its progenitor supernova, and the surrounding environment. In 
this paper, we present a new semi-analytic model for the evolution of such a pulsar 
wind nebula. This model couples the dynamical and radiative evolution of the pulsar 
wind nebulae, traces the evolution of the pulsar wind nebulae throughout the lifetime 
of the supernova remnant produced by the progenitor explosion, and predicts both 
the dynamical (e.g. radius and expansion velocity) and radiative (radio to TeV 7-ray 
spectrum) properties of the pulsar wind nebula during this period. As a result, it is 
uniquely qualified for using the observed properties of a pulsar wind nebula in order 
to constrain the physical characteristics of the neutron star, pulsar wind, progenitor 
supernova, and surrounding interstellar medium. We also discuss the expected evolution 
for a particular set of these parameters, and show that it reproduced the large spectral 
break observed in radio and X-ray observations of many young pulsar wind nebulae, and 
the low break frequency, low radio luminosity and high TeV 7-ray luminosity, and high 
magnetization observed for several older pulsar wind nebulae. The predicted spectrum 
of this pulsar wind nebula also contains spectral features during different phases of its 
evolution detectable with new radio and 7-ray observing facilities such as the Extended 
Very Large Array and the Fermi Gamma-ray Space Telescope. Finally, this model has 
implications for determining if pulsar wind nebulae can inject a sufficient number of 
energetic electrons and positrons into the surrounding interstellar medium to explain 
the recent measurements of the cosmic ray positron fraction by PAMELA and the cosmic 
ray lepton spectrum by ATIC and HESS. 
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Introduction 



The gravitational collapse of the core of a massive star into a neut ron star (e.g. 



Baade &: Zwickv 



19341 ) releases enough energy to power a supernova explosion (e.g. IZwickvlll938l ). The resultant 
neutron star is born spinning, and thought to lose its rotational energy t hrough a ultra-relativistic 
magnetic and particle out flow commonly referred to as a pulsar wind (e.g. iGoldreich &: Julianlll969l : 
Kennel k Coronitil[l984al ^ . The interaction between the pulsar wind and neutron star's environment 
creates an object called a pulsar wind nebula (PWN). Initially (< 10^ years after the supernova 
explosion), the neutron star and its PWN are inside the supernova remnant (SNR) created by the 
expa nsion of the supernova ejecta into the surrounding interstellar medium (ISM). Previous work 
(e.g. [Pacini &: Salvatilll973l : lRevnolds &: Chevalieiill984l ) has demonstrated that, during this period, 
the properties of the PWN depends on the physical properties of the central neutron star (e.g. its 
space velocity, initial period, surface magnetic field strength, spin-down timescale, and braking 
index), the pulsar wind (e.g. its magnetization and energy spectrum of particles injected into the 
PWN), the progenitor supernova explosion (e.g. the mass and initial kinetic energy of the ejecta), 
and the surrounding ISM (e.g. the density profile). Measuring these quantities is important for 
understanding the underlying physical mechanism behind core-collapse supernova, the formation 
of neutron stars in these explosions, and the evolution and properties of the progenitor star. Many 
of these quantities are extremely difficult to measuring directly for most neutron stars / SNRs, but 
possible to infer indirectly using observations of PWNe. 

Using the observed properties of a PWN inside a SNR to determine the properties of the 
neutron star, pulsar wind, proge nitor supernova, and ISM requires understanding its evolution. As 
summarized in a recent review by lGaensler &: Sland (|2006l ). this is extremel y complicated due to the 



rapid evolution of both the SNR and central neutr on star. Analytical (e.g. lOstriker &: Gunn 



Pacini fc Salvatll973l:lE,evnolds fc Ghevaliei] ll984l') and numerical simulations (e.g. 



200 



i 



van der Swaluw et a 



1971 



Bucciantini et al 



20041') of the evolut i on of a PWN inside a SN R identify three important 



Gelfand et al 



2007 



Gaensler &: Sland l2006l ) : an initial free-expansion. 



evolutionary phases (e.g. 
the eventual collision between the PWN and SNR reverse shock which causes the PWN to contract 
and re-expand inside the SNR - eventually stripping the neutron star of its PWN, and a two PWNe 
phase inside the SNR - a "relic" PWN containing the particles injected into the neutron star at 
earlier times and a "new" PWN composed of particles injected by the neutron star after it leaves 
the "relic" PWN. Not surprisingly, the observable properties of a PWN vary significantly during 
this evolution, and are highly sensitive to the physical characteristics of the central neutron star, 
progenitor supernova, pulsar wind, and ISM listed above. Therefore, in order to use the observed 
properties of a PWN to constrain these quantities, it is necessary to have a model for the evolution 
of a PWN inside a SNR which takes into account the relevant physics of all of these components. 



Many such models exist in literature (e 


g. Blondin et a 


2001: 


Bucciantini et al. 


2003, 


2004 




Chevaher & Pransson 199^: 


Chevahei 2005: 


Del Zanna et al. 


2004: 


Gelfand et al. 


2007 


: Jun 


1998 


Kennel & Coronitll984a b: 


Ostri 


cer & Gunn 


1971: Pacini & Salvati 


1973: Eees & GunJl974 


: Revnolds &: Chevali 


19841: van der Swaluw et al. 


2001: 


van der Swaluwiboosi: van der Swaluw et al. 2004 


: Venter & de Jager 
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20061 : IVolpi et al.ll2008l and references therein). While most reproduce the general evolutionary se- 
quence outlined above, they differ significantly in the details, e.g. their treatment of the evolution 
of the PWN's magnetic field, the injection of the pulsar wind into the PWN. These differences 
have significant consequences for the predicted evolution of the observable properties of a PWN 
(such as its size and broadband spectrum). A dditionally, most models pred ict either the dynamical 



properties (e.g. size and e xpansion velocity ; Ivan der Swaluw Sz Wul I2OOII ) or spectral properties 



(e.g. broadband spectrum: IVolpi et al.ll2008l ) of the PWN but not both. In this paper, we present 
a new semi-analytic model for the evolution of a PWN which predicts both the dynamical and ra- 
diative properties of a PWN inside a non-radiative SNR through the entire evolutionary sequence 
described above. This is important because there are many syst ems of interest whe re the PWN 
likely has col l ided with the SNR reverse shock (e .g. G328.4+0.2. ICelfand et ahlboOTl : G327.1-1.1, 
Temim et al.l |2009j; Vela X iLaMassa et al.l |2008| ) . This model also self-consistently couples the 
dynamical and radiative evolution of the PWN, including the evolution of the PWN's magnetic 
field, and calculates the broadband (radio - TeV 7-ray evolution) of the PWN. As a result, it is 
is well suited for both examining the effect of different supernova, neutron star, pulsar wind, and 
ISM properties on the evolution of the resulting PWN, and for using the observed properties of a 
PWN to determine the physical properties of the progenitor supernova, central neutron star, and 
surrounding ISM. 

This paper is structured as follows. In ^ we describe the physics underlying our model and its 
implementation. In ^ we present and discuss the predicted evolution of a PWN for a particular 
set of input parameters. Finally, in ^ we discuss the implication of these result and potential 
applications of this model. 



2. Model Description and Implementation 



In this Section, we describe the underlying physics of this model for the evolution of a PWN 
inside a SNR ( ^2.ip and its implementation ( ^2.2p . 



2.1. Model Description 



This model assumes the PWN is surro unded by a thin shell of material with radius Rpwn, niass 



sw,pwn 



and Gxpaiision velocity I'pwn 

(e.g. lOstriker &: GunnI 



Gelfand et al.l 120071 ). If Vp„n is 



larger than the velocity of the material surrounding this PWN, Ve]{Rpwn), we assume that Mgw.pwn 
increases by an amount AMgw.pwn = 47ri?p^^AiipwnPcj(-Rpwn)) where Pe){Rpwn) is the density of 
material surrounding the PWN. The difference in pressure between the PWN interior to this mas s 



shell, Ppwn, and surrounding SNR, i-*snr(-Rpwn)) applies a force, Fap, equal to (jGelfand et al.ll2007l ): 



FaP — 47ri2p^n [Ppwn — -Psnr(-Rpwn)]- 



(1) 
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The resultant change in momentum of this shell of material is simply: 



di 



(M 



sw,pwn '^pwn J 



Fap, 



(2) 



and we use this equation to determine the dynamical evolution (e.g. i?pwn(i) and Vp^^{t)) of the 
PWN. This approach requires modeling the SNR's density, velocity, and pressure profile with time 
and the pressure inside the PWN, Ppwn- We evolve the properties of the SNR using the results of 
previously developed analytic models described in Appendix [Xj Since no such descriptions exist 
for a SNR in the radiative phase of its evolution, our model only predicts the evolution of a PWN 
inside a SNR during its free-expansion and Sedov- Taylor evolutionary phases. 

The evolution of Ppwn depends on the injection rate of energy into the PWN by the central 
neutron star, the content of the pulsar wind, and the evolution of the particle and magnetic com- 
ponents of the PWN. We assume that all of the spin-down lumi nosity E of the neutron star is 
injected into the PWN (e.g. Equation 5 in lGaensler &: Slandl2006l ): 



m = Eo 



i + — 

Tsd 



P+1 
p-1 



(3) 



where p is the the braking inde^l^, £"0 is the initial spin-down luminosity of the neutron star, and 
Tsd is the spin-down timescale. As mentioned in 511 it is believed that neutron star will eventually 
be stripped of the PWN. When this o ccurs, it no longer injects energy into the "relic" PWN 
but forms a new PWN insid e the SNR ( van der Swaluw et al.ll2004l ). Based on the simulations of 
van der Swaluw et al.l (|2004l ). we assume the neutron star leaves the PWN only after the PWN has 
collided with the SNR reverse shock, when the distance the neutron star has traveled since the 



supernova explosion. 



psr; 



satisfies rpsr > i?, 



"-psr' 



The energy of the pulsar wind is distributed between the electrons and positrons (£'inj,e)5 
(E'inj.i), and magnetic fields (£'inj,B) that comprise this outflow, such that: 



E 



-^inj,e ~l~ E'm^ i + -E'inj,B' 



To parameterize the content of the pulsar wind, we define the following variables: 

Em},c{t) 



Ve{t) = 
Vi{t) = 



m 

Eini,i{t) 

m 

EinlB{t) 



10ns 

(4) 

(5) 
(6) 
(7) 



^The braking index p is defined as tlpsi = ^ 

where P is the rotational period of the neutron star; e.g. Shapiro fc Teukotskv 19861 ) 



fcfipsr wher e Hpsr is the spin frequency o f the neutron star {Q.ps 



2'k/P 
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where r]f> + rji + ijb = ^ ■ It is important to emphasize the r/B is not equivalent to the inagneti zation 



parameter, a, defined as the ratio of magnetic to particle energy (jKennel &: Coronitilll984al ). The 



the ram pressure of the unshocked wind is equal to Ppwn (e 
occur s at a distance from the neutron star, rts, equal to (e.g. 

2nnd ): 





Goldreich &; Julian 


1969 


). This 


Slane et al. 


2004 


1 


Gaensler & Slane 



E 



where ^ is the filling factor of the pulsar wind (an isotropic wind has ^ = 1). In this model, we only 
concern ourselves with the content of the pulsar wind just downstream of the termination shock 
(r > rts), which is likely very different than the content of the pulsar wind near the neutron star 



(e.g. lAronsll2007l ) 



Since the energy of the pulsar wind is divided among electrons and positrons, ions, and mag- 
netic fields, the same must be true for the energy inside the PWN, Ep^^: 



E, 



pwn 



(9) 

is the kinetic 

energy of ions in the PWN, and -Epwn.B is the energy stored in the magnetic field of the PWN. 
Each component contributes separately to the total pressure inside the PWN, such that -Ppwn is: 



where -Epwn e is the kinetic energy of electron and positrons in the PWN, E, 



pwn,i 



Pr 



pwn 



-fpwn,e ~l~ -fpwiiji ~l~ -fpwn,B) 



(10) 



where -Ppwn,e is the pressure associated with electrons and positrons in the PWN, Ppwn,i is the 
pressure associated with ions in the PWN, and Ppwn.B is the pressure associated with the PWN's 
magnetic field. The energy and associated pressure of these components evolve differently from 
each other, as explained below. 

In this model, we assume the magnetic field inside the PWN is uniform and isotropic. As a 
result, -Epwn.B is: 



E, 



pwn,B 



d2 



pwn 



(11) 



where -Bpwn is the strength of the PWN's magnetic field, and the pressure associated with this 
magnetic field Ppwn,B is: 



pwn,B 



d2 

^pwn 

Stt 



(12) 



We also assume that the magnetic flux of the PWN is conserved as i?pwn changes. As a result, 
-Bpwn oc i?p^n neglecting any input from the neutron star. In this case, ^'pwn.B oc iip^n from 
Equation [TT] and -Ppwn,B oc R^^^ from Equation [T21 

We additionally assume that the electrons, positrons, and ions inside the PWN are relativistic. 
Since the density inside the PWN, it is same to assume that these particles are essentially colli- 
sionless and therefore behave as an ideal gas with adiabatic index 7 = 4/3. Therefore, the pressure 
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associated with particles inside the PWN Ppwn,p is: 



pwn,p 



pwn,p 



pwn 



(13) 



where -Epwn.p is the total particle energy of the PWN (-Epwn,p = -£'pwn,e + -E'pwn.i) and l^wn is the 
volume of the PWN (Fpwn = ^T^Rp^n)- The evolution of -Epwn.p and -Ppwn.p depends on adiabatic 
and radiative losses suffered by the particles inside the nebula. Since we are assuming that all 
the particles inside the PWN are relativistic (all three populations has 7 = 4/3), the adiabatic 

evolution of -Epwn.p is independent of the distribution of particle energy between these species as 

4/3 

well as their energy spectrum. If the evolution of the PWN is purely adiabatic, then -Ppwn.p^pwn is 
constant, implying that £'p„n,p oc R~^^ and Ppwn.p oc R~^^. 

The total radiative losses of electrons/positrons and ions in the PWN depends greatly on the 
energy distribution and is sensitive to their individual energy spectra. Their energy spectrum is 
highly sensi tive to the injectio n spectrum of par ticles into the PWN. Obs ervations of several PWNe 
(e.g. 3C58; Isiane et allbood . PWN GO. 9+0.1; IVenter fc de Jaeerl bood l suggest a br oken power- 
law i njection spectrum for electrons and positrons, and recent theoretical work (e.g. ISpitkovskv 
20081 ) argues the injection spectrum of these particles is a Maxwellian with a non-thermal tail. For 
simplicity, we assume a simple power law injection spectrum for the electrons, positrons, and ions, 
but our formalism can easily accommodate a more complex input spectrum. 



Assuming that the injection spectrum of electrons and positrons obeys a single power-law, we 



have: 



E \ 

ne = no^e [ J electrons s^^ keV^^, 



(14) 



where UeAEAt is the number of electrons and positrons with energy between E and E + AE 
injected into the PWN in time At. For consistency with Equation HI we require that: 



E, 



mj,c 



Uf-EdE, 



(15) 



where -Eejinin and -Ee.max are, respectively, the minimum and maximum energy of electrons and 
positrons injected into the PWN. Therefore, the rate of electrons and ions injected into the PWN, 
A^inj,e, is: 



mj,e 



-7e 



-7e 



-^e,max 



^e,min T7' 



HE, 
HE. 



e,max 
e,max 



e,m 
e.mm 



E- ■ 



E. 



7e 
7e 



(16) 



-'^e.max -'^e.rnin 

Similarly, we assume the injection spectrum of ions is a single power law: 



no,i 



E 
Eq 



(17) 
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where rijAi^At is the number of ions with energy between E and E + AE injected into the PWN 
in time At. Again, for consistency with Equation [J] we require that: 



E, 



TiiEdE, 



(18) 



where -Ei^min and Ei^^g^^ are the minimum and maximum energy of ions injected into the PWN. 
This results in the injection rate of ions into the PWN, Mnj^j, having the same functional form as 
Equation [16l 

This model assumes the dominant radiative processes are synchrotron emission and inverse 
Compton scattering off background photons. For a particle with energy E, mass m, ar id charge q. 



the r ate at which it loses energy due to synchrotron emission, -Psynch) is (Equation 3.32 in lPacholczvk 
197d ): 

Psynci. = ^^Bl^^sm'eE\ (19) 

where 9 is the angle between the particle's velocity and the magnetic field. We assume that the 
velocities of p articles in the PWN are raii domly oriented with respect to the magnetic field, such that 
sin^ ^ = 2/3 ( Rybicki Lightman 197S). W hen the PWN is small, synchrotron self-absorption is 
important (e.g. [Reynolds &: Chevalieiill984l ). To determine its effect on the dynamics and emission 
of the PWN, we first assume that a particle emit s all of its synchr otron radiation at a frequency 
u = Ucrit, where Vciit is equal to (Equation 3.28 in iPacholczvkl 119701 ): 



^crit 



-. — sin ( 



(20) 



where sin 9 = -y/2/3. The optical depth of the PWN to this emission is (§3.4 in IPacholczvkl Il97d ) : 

-ppRp^^K{l), (21) 



, ^ ,0 ,.7(^PwnSin0)3/2 

-r(i^crit) = 4.2 X 10' — ^ 



5/2 
^crit 

where sin^ = -\/ 2/3 as before, pp is the density of particles inside the PWN, and K{1) = 1.1 



(|Pacholczvklll970l ). We assume that particles with energy such that r(t/crit.) > 1 do not lose e nergy 
through synchrotron radiation - an approach similar to that of lRevnolds &: Chevalieii (| 19841 ). 



The power lost by particles due to inverse Compton scattering off background photons, Pic, 
is (Equation 5.57 in |Pacholczvklll970l ) : 



Pic 



9(mc2)4 



u,^dE^f{E), 



(22) 



where Urad is the energy density of the background photon field, and f{E) is the inverse Compton 
scattering cross-section relative to the Thomson cross-section. f{E) is sensitive to both the particle 
energy and the spectrum of the background photon field. The background photons are believed to 
be dominated by the Cosmic Microwave Background (CMB; Uj-ad ~ 

4.17 X 10 



13 gj.gg 3^^ g^j^^ 
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we calculate f{E) for this field using the procedure described in §2.3 of lVolpi et al.l (|2008l ) (Figure 
[1]). It is possible that other photon fields, e.g. starlight, not considered here are important for a 
given PWN - and a simple adjustment to the formalism presented here can account for these as 
well. 



2.2. Model Implementation 

To determine the evolution of a PWN inside a SNR, we use the equations presented in ^2.11 to 
determine how the relevant properties of a PWN at a time t evolve to a time t + At. The procedure 
used to determine the initial conditions of the PWN inside a SNR is described in Appendix [Bl and 
the free parameters of this model are listed in Table [TJ We calculate the properties of the PWN at 
time t + At using the following procedure: 

• Step 1: Calculate the new radius of the PWN 

Rp^nit + At) = i?pwn(t) + t'pwn(t) X At, (23) 

and the distance the neutron star has traveled inside the SNR since birth: 

rpsr(t + At) = Vpsr{t + At) X (t + At), (24) 
where Vpsi is the space velocity of the neutron star. 

• Step 2: Adjust the electron/positron, ion, and magnetic field energy of the PWN for the 
adiabatic work done by the PWN on its surroundings and, if the neutron star is still inside 
the PWN, the energy it injected between times t and t + At: 

R (t) /-i+At _ 

£;pwn,e(t + At) = „ 7, \,, Ep^nAt)+r]e Edt (25) 

Ep^r.,i{t + At) = Ep^^,i{t)+r^i r^^' Edt (26) 

^pwnit + Atj Jt 



Rpwn(i} 
Rp^n{t + Aty 



Spwn,B(t + At) = ^^^ii^Ep^„^B(t)+r/B / Edt (27) 



where E is defined in Equation [3] and J Edt is solved analytically. We then use £'pwn,B(i + At) 
to calculate -Bpwn(i + At) using Equation [11] and Ppwn,B{t + At) using Equation [T2l 

• Step 3: Calculate the energy spectrum of electrons, positrons, and ions inside the PWN 
at time t + At taking into account adiabatic, synchrotron, and inverse Compton losses. A 
particle with energy E at time t has an energy at time t + At equal to: 

+ = RpUt^+^tf " ^^^""^'^^^ ^ ^c(t)]At (28) 
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where -Pgynch (t) (Equation [T9|) is calculated using i3pwn (t) and Pic (t) is defined in Equation 
[22l with the restriction that E{t + At) > 0. We then subtract from the value of Ep^^^eit + At) 
and Ep^^^i{t + At) calculated in Step 3 the total synchrotron and inverse Compton losses of 
these particles. Even though the radiative losses of ions in the PWN is significantly smaller 
than that of electrons and positrons (both Pgynch and Pic are oc ), for completeness they 
are calculated. We then determine Ppwn,pit + At) using Equation [T3l 

• Step 4: Calculate the properties of the surrounding SNR [pej{Rp^^,t + At), VejiRpwm t + At) , 
and Psnr(-Rsnr, t + At)] using the prescription described in Appendix lAl If Vej{Rpwn, t + At) < 
fpwn(t) then: 

Msw,pwn(t + At) = Msw,pwn(t) + ^VT [i?p„n(t + At)^ - i?p„n(t)^] Pcj(i?pwn, t + At). (29) 

Otherwise, Msw,pwn(i + At) = Ms„,pwn(i)- 

• Step 5: Using the properties of the PWN and SNR determined above, calculate F/\p{t) 
(Equation [1]) . The new velocity of the mass shell surrounding the PWN, t'pwn(i + At) is : 

Vpy,n{t + At) = — ^ , AA ^ [^sw,pwn(^)^'pwn(^) + AMsw,pwn^^ej(-Rpwn,^)+i^pwn(^) X At], 

-'"sw,pwn(,t + ZAtj 

(30) 

from conservation of momentum arguments, where AMsw,pwn = Afsw,pwn(i+ At) — Msw,pwn(i)- 



3. Model Performance 

In this Section, we present the evolution of a PWN inside a SNR predicted by this model 
for the combination of the input variables listed in Table [2l It is important to emphasize that 
these results are specific to this set of parameters, and the different PWNe may have very different 
values for some or all of these quantities. The properties of the progenitor supernova {Eg^ and 
Mej), central neutr o n sta r (£"0, Tgd, and p), and surrounding ISM (nism) are the same as Model 



A in iBlondin et al.l (|200ll ) - chosen by these authors to reproduce a PWN similar to that of the 
Crab Nebula. We also assume a purely electron/positron pulsar wind (i.e. r/j = 0), and that the 
properties of the pulsar wind (r/^, rj^., i?e,min) -E'e,max) 7e) remain constant with time. We also use a 
pulsar velocity of 120 km s ~^, the most recent measurement of the transverse velocity of the Crab 



pulsar (jKaplan et al.l 120081 ) . As shown in Figure [2l for this particular set of input parameters our 



model predicts four evolutionary stages: 

• Initial expansion: This phase ends tcoi ~ 4500 years after the supernova explosion when 
the PWN collides with the reverse shock. 

• Reverse Shock Collision and First Compression: This phase ends ~ 20000 years after 
the supernova explosion. During this phase, ~ 17000 years after the supernova explosion, we 
expect the PWN to be stripped of its neutron star. 
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• Re-expansion: This phase ends ~ 56000 years after the supernova. During this phase, the 
neutron star re-enters the "rehc" PWN created during the first compression ~ 30000 years 
after the supernova explosion. 

• Second Compression: This phase continues until the SNR enters the radiative phase of its 
evolution. During this contraction, the PWN is stripped of its neutron star ~ 70000 years 
after the supernova explosion. 



The dynamical and radiative evolution of the PWN in these four phases is described in detail below. 

The initial expansion of the PWN is the result of Ppwn ^ -fsnr(-Rpwn) (Figure [3]). This is the 
case because Psnr(-Rpwn) ~ when .Rpwn < Rrs, i-e. adiabatic expansion of the SNR makes the 
ejecta interior of the reverse shock cold (Appendix |A|) . During this time, the PWN is expanding 
faster than the surrounding ejecta (Figured]), and Msw.pwn increases (Figure [5]). For most of this 
expansion, the energy input for the neutron star (E) is larger than the adiabatic losses resulting 
from the PWN's expansion and the radiative losses suffered by electrons inside the PWN (Figure 
E]), causing the internal energy of the PWN to rise (Figure [Tj). At the very earliest times (t < 100 
years), synchrotron losses approach the energy input of the pulsar due to the PWN's very strong 
magnetic field (Figure[8]). These synchrotron losses cause the magnetization parameter of the PWN, 
a, to temporarily increase by more than an order of magnitude (Figure [9|). For t < 1000 years, 
synchrotron losses actually dominate over adiabatic and inverse Compton losses (Figure [6]), causing 
the PWN to expand slower than predicted by Equation IB7I This difference is small at early times, 
so Equation IB7I does provides a reasonably accurate initial condition. The expansion of the PWN 
causes -Bpwn to decrease rapidly (Figured]), resulting in a rapid decrease of the PWN's synchrotron 
luminosity. As a result, adiabatic losses dominate from ~ 1000 years after the supernova explosion 
until the PWN collides with the reverse shock. These adiabatic losses cause both -Bpwn and -Epwn,B 
to decrease for t > 1000 years, while -Ep„n^p increases due to the continued inject of particles into 
the PWN by the pulsar. It is im portant to note that -Bp„ ri does noi follow .Bpwn oc 1/(1 + (t/rsd)"), 



as assumed in other work (e.g. IVenter &: de Jageii 120061 ). We find that -Bpwn cx: t ^ durin g the 



initial expansion, similar to the t — t ^ behavior derived bv lRevnolds Chevalieii (jl984l ) 



This behavior ends when the PWN collides with the SNR reverse shock. This collision shocks 
the swept-up ejecta surrounding the PWN, but not the PWN itself because the sound speed inside 
the PWN (~ c/Vs) is significantly higher than the velocity of the reverse shock. The collision 
with the reverse shock marks the end of PWN's rapid expansion (Figured]) because the PWN is no 
longer in an essentially pressureless environment. At the time of this collision, Ppwn -fsnr(-Rpwn) 
(Figure [3]), causing Vp^^ to decrease significantly (Figure H]). However, for the first ~ 100 years 
after this collision, Vp^n > i'ej(-Rpwn) (Figure H]), resulting in Msw,pwn continuing to rise, increasing 
from ~ IMq to ~ 3Mq (Figure [5]) due to the high density behind the reverse shock. The sharp 
decrease in Vpwn also leads to a sharp decrease in adiabatic losses (Figure [6]) . Because the pulsar 
continues to inject energy into the PWN, both E'pwn.p and -Epwn.B initially rise rapidly after the 
PWN/reverse shock collision (Figure [7]), though Ppwn is still Ppwn ^ ^snr(-Rpwn)- As a result, the 
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high-pressure ejecta downstream of the reverse shock will compress the PWN beginning ~ 6500 
years after the supernova explosion. 

This compression begins a series of contractions and re-expansions that continue until the SNR 
enter s the radiative phase of its evolution (Figure [ 2l). Similar behavior was observed in previous work 
(e.g. Blondin et al.ll200ll : iBucciantini et al.ll2004l : Ivan der Swaluw et al.ll200ll . 12004 ). The adiabatic 
compression of the PWN causes a sharp rise in the internal pressure (Figure [3]) and the particle and 
magnetic energy (Figure [7|) of the PWN. The compression of the PWN also causes a rapid rise in 
Bpwn (Figure [8|), which in tu rn results in a rapid ris e in th e PWN's synchrotron luminosity (Figure 



[6]) - as previously noted by iRevnolds &: Chevalieii (|1984| ) . Eventually, the synchrotron losses are 
larger than the work done by the surrounding SNR on the PWN and the rate particle energy is 
injected into the PWN by the neutron star, causing E'pwn.p to decrease (This would have occurred 
during the second compression if we evolved the PWN further in time). This decrease in Ep^^^p 
leads to a decrease in the synchrotron luminosity of the PWN (Figure [6]) , even though the magnetic 
field inside the PWN is still increasing (Figure [8]). Since synchrotron losses do not diminish -E'pwn,B; 
a increases significantly during the compression phases of the PWN (Figure [9]) . 

Due to the space velocity of the neutron star, we expect that the PWN will be stripped of its 
pulsar during both compressions. It is important to note that this result is strongly dependent on 
the space velocity of the pulsar. For a lower space velocity, it is possible that the PWN will not be 
stripped of its pulsar until the pulsar leaves the SNR, while for a higher velocity the PWN will not 
overtake the pulsar during its re-expansion. The departure of the pulsar causes a rapid increase of 
a (Figure [9]) and a rapid decrease in the synchrotron luminosity of the PWN since it removes the 
only source of high energy particles that dominate this emission (Figure [6]) . The latter does not 
occur during the second compression due to the rapid rise in Spwn when the pulsar exits the PWN. 
It is expected that the pulsar cre ates a new PWN inside the SNR after it leaves its original PWN 



(e.g. Ivan der Swaluw et al.ll2004l ) whose properties we do not calculate. 



Eventually during the compression stage, the pressure inside the PWN will exceed that of 
the surrounding SNR, causing the PWN to re-expand. During the first contraction, the pressure 
inside the PWN exceeds that of the surrounding SNR ~ 5000 years after the compression begins 
(Figure [3]). At this time, the PWN is contracting at a very high velocity (> 1000 km s~^; Figure 
H]) and, due to the considerable momentum of the swept-up material surrounding the PWN, it 
takes ^ 8000 years for PWN to re-expand (Figure [2|). Therefore, Ppwn ^ ^snr(-Rpwn) when re- 
expansion begins (Figure [3|). During this phase, the PWN sweeps up additional material. Due to 
the low density inside the SNR at this time, the mass shell only accumulates ~ 0.5Mq of additional 
material. Initially during the re-expansion, -Epwn.p) -E'pwn.B (Figure [6|), and .Bpwn (Figure [8|) rapidly 
decrease due to adiabatic losses. The decrease in -Epwn and increase in -Rpwn result in a rapid 
decrease in -Ppwn (Figure [3|), and the decrease in both -Bpwn and -Epwn.p result in a rapid decrease 
of the PWN's synchrotron luminosity (Figured]). During the re-expansion, a initially continues to 
increase because the synchrotron and adiabatic losses suffered by particles in the PWN are greater 
than the adiabatic losses of the magnetic field (Figure [9]) . 
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During the re-expansion, Wpwn reaches ~ 300 km s~^, higher than the space velocity of the 
neutron star. As a result, the pulsar will re-enter the "relic" PWN. At the time of re-entry E is 
actually higher than the synchrotron, adiabatic, and inverse Compton losses of the PWN (Figure 
[6]) , causing -Epwii to increase (Figure [7]) . Since most of the energy deposited by the pulsar in the 
PWN is in the form of particles (r/g r/B ) , only ii^pwn.p increases as a result of the re-entry of the 
pulsar - -Einj.B is significantly less than the adiabatic losses suffered by the magnetic field (Figure 
[7]). As a result, a begins to decrease, though cr » 77^ at these times. The increase in Spwn.p also 
leads to a sharp rise in both the synchrotron and inverse Compton luminosity of the PWN (Figure 
El . This rise in energy is not accompanied by a rise in pressure (Figure ED because the volume of the 
PWN is increasing faster than the particle energy. As a result, ~ 30000 years after the supernova 
explosion Ppwn < -fsnr(-Rpwn), and the PWN will be compressed for a second time. As shown in 
Figure [21 this second compression is expected to continue until the SNR enters the radiative phase 
of its evolution, and the evolution of the PWN during this phase is described above. 

The complicated dynamical evolution of a PWN described above causes the energy spectrum 
of electrons and positrons in the PWN, and the photons they radiate, to change considerably 
with time. As mentioned earlier, synchrotron losses play a very important role during the initial 
expansion of the PWN. For t <C r^d, the magnetic field of the PWN is so strong that the only high 
energy {E > 1 TeV) electrons and positrons in the PWN are those which were injected very recently 
(Figure [TOj) . For r^d < t < tcoh the electron and positron spectrum of the PWN is dominated by 
previously injected particles at all energies, and the sharp high-energy cutoff observed at early 
times is replaced with a more gradual turnover whose sharpness increases with time (Figure [TO]l . 
This is due to the increase peak energy of the electron and positron spectrum (as first noted by 
Reynolds &: Chevalieiill984l ) resulting from the decrease in -Bpwn- 



Correspondingly, the spectrum of photons radiated by the PWN during its initial expansion 
also evolves with time (Figure [TT]) . As explained in ^ the only emission mechanisms we consider 
are synchrotron radiation and inverse Compton scattering off the CMB, and the photon spectrum 



is calculated using the procedure described in Section 2 of lVolpi et al.l (120081 ). At very early times 
(t < 50 years), the PWN is extremely luminous at GeV energies due to the large number of high 
energy electrons and positrons recently injected by the neutron star (Figure [T2]) . For t < Tsd, the 
PWN's synchrotron emission peaks at photon energies > 100 keV (Figure [TT]) . and the PWN's 
luminosity is highest in the hard X-ray regime (Figure [T2]) . Since recently injected particles are 
energetically important, both the synchrotron and inverse Compton spectrum have two peaks - a 
lower energy peak resulting from previously injected particles and a higher energy peak resulting 
from recently injected ones. For Tgd < t < t^oi, the synchrotron and inverse Compton spectrum have 
a single peak since recently injected particles no longer dominate. During the initial expansion, 
the energy peak of the synchrotron spectrum decreases from ~ 1 MeV at t <C Tgd to ~ 1 — 10 keV 
at t ~ tcoi, causing a rapid decrease in the hard X-ray luminosity of the PWN. Slower decreases 
are predicted for the radio - soft X-ray luminosity (Figure [T2]) due to the gradual decline in the 
synchrotron luminosity of the PWN (Figure ^ resulting from the decreasing value of .Bpwn (Figure 



-IS- 



IS]). Conversely, the energy peak of inverse Compton emission increases from ~ 1 TeV at i ~ 
to ~ 50 TeV at t ~ tcoi due to the increase in the break energy of the electron spectrum discussed 
above, leading to a rapid rise in the GeV and TeV 7-ray luminosity of the PWN (Figure [T2]) . 

The energy spectrum of electrons and positrons inside the PWN changes significantly during 
the first contraction (Figure [TO]) . Due to the strong magnetic field inside the PWN during this phase 
(Figure [8]), the synchrotron lifetime of the highest energy (E>10 TeV) electrons and positrons in 
the PWN becomes significantly less than the age of the PWN. This results in a sharp cut-off in the 
electron and positron energy spectrum where the synchrotron lifetime of electrons and positrons 
is the age of the PWN, above which recently injected particles dominate. The strengthening -Bpwn 
results in this energy decreasing with time. This causes the peak photon energy of the synchrotron 
emission to decrease significantly during this time, from ~ 10 keV when the PWN collides with 
the reverse shock to ~ 100 eV when the neutron star exist the PWN. The spectrum of the inverse 
Compton emission radiated by the PWN also changes considerably during this time - with the 
energy peak decreasing from ~ 10 TeV at the time of the PWN/reverse shock collision to ^ 100 
GeV when the neutron star leaves. The luminosity of the PWN in the wavebands dominated by 
synchrotron emission (radio - soft X-rays; Figure [TT]) increase (Figure [T2]) due to the strengthening 
magnetic field. The hard X-ray and GeV 7-ray luminosity of the PWN increases due to the 
decreasing energy of the inverse Compton peak, which causes a decrease in the TeV 7-ray luminosity 
of the PWN (Figure [T2]) mitigated by inverse Compton emission from the highest energy recently 
injected particles. 

The departure of the neutron star from the PWN during the first contraction has a dramatic 
effect on the electron and positron energy spectrum in the PWN. When this occurs, the only high 
energy {E > 10 TeV) electrons and positrons in the PWN are those recently injected by the neutron 
star. Therefore, the departure of the neutron star causes this plateau of high-energy electrons and 
positrons to quickly disappear due to synchrotron cooling (Figure fTO]) . This results in a sharp 
decrease in the PWN's X-ray and TeV luminosity (Figure [T2|) . It is possible that the new PWN 
will be luminous at these energies, so the total X-ray and TeV luminosity of the system might be 
appreciably higher than predicted here. The energy peak of the electron and positron spectrum 
also decreases from ~ 10 TeV when the neutron star leaves to ~ 10 GeV at the time of 
re-expansion (Figure [TO]) . with the energy spectrum highly suppressed above the peak. This causes 
the peak of the synchrotron emission to decrease from the ultra-violet to mid-IR wavelengths and 
the bulk of the PWN's synchrotron emission to move from the optical to radio wavelengths, while 
the peak of the inverse Compton emission to decrease from ~ 100 GeV to ~ 10 MeV (Figure [TT]) . 
causing the PWN to be more luminous at GeV energies then TeV energies (Figure [T2]) . 

As the PWN re-expands into the SNR but before it overtakes the pulsar, the maximum energy 
of electrons and positrons in the PWN decreases from ~ 10 GeV to ~ 1 GeV (Figure [TO]) . This 
fact, and the decrease in -Bpwn (Figure [8]), causes the peak frequency of the PWN's synchrotron 
spectrum to decrease to ~ 1 GHz and the peak of the inverse Compton radiation to decrease to 
~ 100 keV. As a result, when the PWN is re-expanding inside the SNR but before it overtakes the 
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neutron star, almost all of its emission is in the radio, soft X-ray, and hard X-ray bands (Figure 
[T2|) . When the neutron star re-enters the PWN, it resumes injecting high-energy electrons and 
positrons into the PWN - resulting in an electron spectrum with two distinct components: a lower 
energy population composed of electrons and particle injected at earlier times and a higher energy 
population composed of recently injected particles (Figure [T0|) . During the PWN's re-expansion, 
the peak energy of the "old" particles continues to decrease due to adiabatic losses while the peak 
energy of the "new" particles decreases due to synchrotron losses (Figure [TO]l . This dichotomy 
also extends to the photon spectrum (Figure [TTIl . The radio emission from the PWN is dominated 
by synchrotron radiation from the old particles, and the radio luminosity of the PWN decreases 
(Figure [T2]) due to their decreasing energy. At all other wavebands, the emission is dominated by 
recently injected particles. Initially, the synchrotron emission from these electrons and positrons 
extends to ~ 100 keV, but due to adiabatic and synchrotron losses this decreases to ~ 10 keV by 
the PWN begins its second contraction (Figure [TT]) . As a result, after an initial increase which 
marks the neutron star's re-entry, both the soft X-ray and hard X-ray luminosity of the PWN will 
decrease - though the decline in the soft X-ray luminosity is significantly slower than the hard 
X-ray luminosity (Figure [T2]) . The inverse Compton emission from the recently injected electrons 
and positrons extends to TeV energies (Figure [TT]) , and the 7-ray luminosity of the PWN increases 
between the re-entry of the pulsar and the second compression (Figure [T2]l . 

During the second contraction of the PWN, the energy spectrum of electrons and positrons 
inside the PWN maintains the two component structure created during the PWN's re-expansion 
(Figure [TOl) . The peak energy of the "old" electron population increases slightly from ~ 1 GeV 
to ~ 5 GeV due to the work done on the PWN by the surrounding ejecta, but the energy peak 
of the "new" electron population decreases from ~ 100 TeV to ~ 1 TeV (Figure [TOl) due to the 
increasing strength of the PWN's magnetic field and the second departure of the pulsar from the 
PWN ~ 15000 years after the second contraction begins. During the second contraction, the photon 
spectrum of the PWN has four distinct regimes - from lowest to highest energy, they are synchrotron 
emission from "old" electrons and positrons, the synchrotron emission from "new" electrons and 
positrons, the inverse Compton emission from "old" particles, and the inverse Compton emission 
from "new" particles (Figure [TT]) . During the second compression, the energy peaks corresponding 
to emission from the old particles increase - the peak frequency of the synchrotron emission rises 
from u 10 MHz to 1/ ~ 1 GHz, while the peak energy of their inverse Compton emission rises 
from ~ 10 keV to ~100 keV. This results in an increase in the radio and hard X-ray luminosity 
of the PWN (Figure [T2|). Conversely, the energy peaks of the spectral features produced by the 
recently injected particles decreases - the synchrotron peak decreases from ~ 10 keV to ~ 1 eV, 
while the inverse Compton peak decreases from ~ 100 TeV to ~ 10 GeV, causing an increase in 
the optical and mid-IR luminosity and a decrease in the 7-ray luminosity of the PWN. 

The observable properties discussed above are calculated using the predictions of the model. 
Unfortunately, for most PWNe the only observational data available is a flux density measurement 
at two or three radio frequencies (typically, 1.4, 4.8, and 8.5 GHz) and a measurement of the 
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X-ray spectrum between E.y ~ 0.5 — 10 keV - though due to interstellar absorption, for several 
PWN it is oiily po ssible to measure the spectrum between = 2 — 10 keV (e.g. 0328.4+0.2; 



Gelfand et al.l 120071 ). This information is then used to estimate the properties of the PWN, for 
example its energetics and strength of its magnetic field. The radio and unabsorbed X-ray spectrum 
are usually fit to a power law, L^, oc u", where a is the spectral index. (The X-ray spectrum is 
often fit using oc , where Ni, is the number density of observed photons and F is the photon 
index, where a = 1 — T.) We derive the spectral index of the PWN in the radio, 0.5 — 10 keV, and 
2 — 10 keV bands throughout its evolution, finding that the radio spectral index is —0.3 except 
for two periods where a <C —1.6 (Figure [TB|) . These two period are during the re-expansion of the 
PWN when the frequency of the lower-energy synchrotron peak decreases from u ^ 8.5 GHz to 
u <1A GHz and during the second contraction of the PWN when the low-energy synchrotron peak 
transitions from u < 1.4 GHz to > 8.5 GHz (Figure [T^ . The photon index F measured between 
0.5-10 keV and 2-10 keV is considerably more variable, fiuctuating between F 1 and F ^ 3. 
As a result, the change in spectral index (Aa) between the radio and the two X-ray bands varies 
significantly with time (Figure [TSl) , with Aa oft en Aa > 0.5, the value expected from standard 
synchrotron theory for a constant magnetic field (jPacholczyld Il970l ) . 



There are several different ways to estimate the internal energy and strength of the PWN's 
magnetic field using the radio and X -ray spectrum of the PWN. In this paper, we evaluate the 
method described by IChevalieii (j2005l ) . This procedure requires the radio and X-ray spectral index 
of the PWN, the inferred break frequency z^b) and luminosity density of the PWN at the break 
frequency, L^^, between these bands, and assumes that -E-pwrijE — (3/4)£'p^n p (minimum energy 
estimate). The values of fb and L^,^ extrapolated from the radio and X-ray spectrum of the PWN 
over-estimates L^^^ due to curvature in the spectrum between these two bands (Figure fill) . This 
fact, coupled with the assumption that a = 3/7 - significantly higher than the actual value of 
a (Figure [9]), causes this method to under-predict -Epwn by a factor of ~ 5 — 10 (Figure [T6]) and 
over-predict .Bpwn by factors of a few at most times (Figure [T7|) . 



The above discussion assumes the PWN is not disrupted by any hydro dynamical instabilities 
as it evolves inside the SNR. However, this is not necessarily the case. The shell of swept-up 
material surrounding the PWN is unstable to Rayleigh- Taylor instabilities when Ppwn > ^snr(^pwn) 
because the low densi ty pulsar wind (ppwn) is accelerating this much higher density shell (psw,pwn) 
(j Chandr asekhaiil 1 96 ll ) . The growth rate of these instabilities depends strongly on the magnetic field 
strength inside the PWN parallel to the boundary between the pulsar wind and swept-up supernova 
ejecta, -Bpwn,||- Numerical simulations of the magne tic field inside PWN su ggest that, at early times. 



the PWN's magnetic field is largely toroidal (e.g. Ivan der Swaluwl 12003), in which case B, 



■Bpwn - though polarized radio observations of older PWNe (e.g. Vela X; 



Milnelll98d . lDodson et al 



20031 ) suggest a strong radial component to their outer magnetic f ield. Assuming B 



the growth rate ujrt{k) of a Rayleigh- Taylor (Kruskal-Schwarzschild; iKruskal Schwarzschildlll954l ) 



B, 



pwn 
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instability with wavenumber k = 2tt/X is (jChandrasekhaiill96ll : iBucciantini et al.ll2004l ): 
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H" Ppwn) 



(31) 



where Opwn is the acceleration of the shell of swept-up material (a 



pwn 



dvpv/n/dt). Hydrodynamic 
simulations of the expansion of the PWN inside a S NR suggest that this shell of swept-up material 
that surrounds the PWN has thickness ~ ^-Rpwn ( van der Swaluw et al.ll200ll ). and we assume 
this is true at all times when calculating Psw.p wn- As a result, the max imum wavenumber /ccrit of a 
Rayleigh- Taylor instability which can grow is (jBucciantini et al.ll2004l ): 
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and the wavenumber of the Rayleigh- Taylor ins tability with the highest growth rate /c, 



( Chandrasekhar 



1961 



(32) 



^crit/2 



Stone &: Gardineij 120071 ) . Rayleigh- Taylor instabilities between the pulsar 



wind and the swept-up ejecta result in "bubbles" of swept-up ejecta entering the PWN. Three- 
dimensional numerical simulations of the growth of Rayleigh- Taylor instabilities in a PWN-like 
scenario (e.g., a strong magnetic field in the light fluid parallel to the interface with the heavy fluid) 
suggest that the penetration of bubbles is relatively unimpeded by the pr esence of a magnetic field , 
though the magnetic field suppresses mixing between these two fluids (jStone &: Gardineij 120071 ) . 
Therefore, to estimate the penetration of swept-up ejecta bubbles into the PWN, we use the results 
of relevant labo ratory experiments. These experiments derived th at the penetration depth of these 
bubbles, /irt, is (jPimonte &: Schneiderlll996l : iDimonte et al.ll2007l ): 



h 



rt 



^Japwn{t)dt 



(33) 



where ab is an experimental derived constant {ob 
the Atwood number of this system: 



0.061; iDimonte fc Schneideiill996l ). and A is 



A = 



Ppwn 



PmSjpwn ~l~ Ppwn 



(34) 



Since p 



ms,pwn 



> Ppwn, for a PWN A ^ 1. 



We find that the PWN is unstable to Rayleigh- Taylor instabilities during the initial expansion 
{t < tco\), and parts of the first contraction, re-expansion {t ~ 15000 — 30000 years), and second 
contraction (t > 85000 years). During the initial expansion, we expect that only Rayleigh- Taylor 
instabilities on the smallest angular scales (< 1°) are suppressed by the nebular magnetic field 
(Figure [THj) , while during the first and second contractions Rayleigh- Taylor instabilities at signif- 
icantly larger angular scale are suppressed (Figure [THj) due to the PWN's strong magnetic field 
(Figure [8|). The growth rate of these instabilities depends significantly on their angular scale, and 
during the initial expansion instabilities with an angular size of ~ 6° grow the fastest (Figures [THl 
fc I19p . Initially during the contractions, the Rayleigh- Taylor instabilities with the highest growth 
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rates have an angular scale ~ 5° (Figure [TH] &: [T9|) . However, as the PWN re-expands, the growth 
rate of instabilities on these angular scales decreases and eventually are suppressed - for example, 
Rayleigh- Taylor instabilities with an angular size of ~ 6° can grow for only ~ 5000 years (Figure 
T9|) . As a result, during this period Rayleigh- Taylor instabilities with large angular scales (> 30°) 
likely experience the most growth (Figure [T9|) . Additionally, the growth-rate of Rayleigh- Taylor 
instabilities at all angular scales is much lower after the PWN has collided with the reverse shock 
than during the initial expansion due to the stronger nebular magnetic field and significantly higher 
density of the mas s shell surrounding P WN (Figure [T9]) . Since the density of these perturbations 



grows as p oc e'^"'' (jChandrasekhaijll96lh . during the initial expansion it is likely that a significant 



fraction of the mass swept-up by the PWN will be in these filaments. 

Rayleigh- Taylor instabilities will cause "bubbles" of swept-up material to penetrate the PWN. 
During the initial expansion of the PWN, the depth of these bubbles is /irt = (0.01 — 0.1)-Rpwn 
(Figure [20]l ■ so ~ 5 — 20% of the volume of the PWN is contained inside this penetration layer. It 
is interesting to n ote that this length scale is similar to the size of the optical filaments surrounding 
the Crab Nebula (iHester et al.lll996l ). After the PWN collides with the reverse shock. Equation [331 
predicts that the depth of this mixing layer decreases significantly, with hj-t ~ when the PWN 
begins to contract (Figure [20]) . This mixing layer will grow again when the pressure inside the 
PWN is higher than that of the surrounding SNR. When the PWN re-expands into the SNR, this 
mixing layer is almost as large as the PWN itself (Figure [20]) . This suggests ta t the PWN might be 



2001 



disrupted at this time, as o bserved in previous hydrodynamical simulations (jBlondin et al. 
van der Swaluw et al.ll2004l ). If so, this will inject a total energy of ~ 6 x 10^^ ergs in the form of 



relativistic particles into the surrounding SNR. It is important to reiterate that the calculation of h^t 
used here (Equation I33p ignores any damping the PWN's magnetic field might have on the growth 
of this maxing layer - as do the hydrodynamical simulations cited above. Recent simulations of 
the growth of Rayleigh- Taylo r instabilities inside a P WN suggest the growth of this mixing layer is 
highly suppressed for high a (jBucciantini et al.ll2004l ). Since this is the case during the contraction 
of the PWN (Figure [9|) , it is likely this approach under-predicts the lifetime of a PWN inside a 
SNR. 



4. Discussion and Conclusions 

In this paper, we present a general model for the evolution of a PWN inside a SNR (f|2|) 
and the specific evolution predicted this model for a particular set of neutron star, pulsar wind, 
supernova explosion, and ISM properties (33]). As mentioned in ^ the ultimate goal of this model 
is to reproduce the observed dynamical and radiative properties of a PWN in order to study the 
central neutron star, progenitor supernova explosion, and pulsar wind. This requires that our model 
accurately reproduces the observed properties of a well-studied and constrained PWN. The best 
test case is the Crab Nebula, the brightest radio and X-ray PWN in the Milky Way and whose age 
(« 950 years old) and neutron star properties {Eq, p, Tgd) are well known. As mentioned earlier, the 



- 18 - 



input parameters of the simulation discussed in ^were based on previous analyses of this source. 
As shown in Table O this model is able to reproduce the size of the Crab Nebula, the radius of the 
termination shock, the expansion velocity of the PWN, and the spectral index of both the radio 
and 0.5-10 keV emission from the PWN. However, for the set of input parameters given in Tabled 
this model predicts a radio. X-ray, and TeV 7-ray luminosity ~ lOx different than that observed. 
Since the relationship between the properties of a PWN at a given time and the input parameters 
of this model is non-trivial, determining if this model can reproduce all of the observed properties 
of the Crab Nebula requires a much more thorough examination of the possible parameter space, 
which we leave for future work. This situation is true for any PWN, and not just the Crab Nebula. 

The evolution of the particular PWN discussed in ^ does have some interesting implications 
for existing questions in the field. The spectrum of this PWN does show the large spectral breaks 
inferred from rad io and X-ray observ ations of several young PWNe though, though not their low 
break frequency (IWoltjer et al.l 119971 ). By exploring the parameter space of possible neutron star, 
pulsar wind, supernova, and ISM properties, we will be able to identify what regimes are necessary 
to produce a low break frequency and large spectral break between the radio and X-ray wavebands 
at early times. For the set of parameters presented in ^ this model predicts a very low break 
frequenc y ^ 1 GHz) at late times - as observed in some older (t > 10^ years) PWNe (e.g. 
DA 495; iKothes et al.l 120081 ). This model also predicts that this particular PWN will have a low 
radio luminosity but a hig h TeV 7-ray luminosity at late times, similar to several PWN recently 
discovered by HESS (e.g. lAharonian et al.l l2006l ) . Additionally, this model that - for this set 
of parameters ~ the magnetization of the PWN a is considerably higher after the PWN collides 
with the reverse shock than before. This could explain why the value of a estimated for the Vela 



PWN (0.05 < cr < 0.5 : ISefako fc de Jager 



LaMassa 



et alJ 



2OO3I ) . a system where this collision is believed to have 



20081 ) . is considerably higher than that of the Crab Nebula 



already occu rred (e.g. 

(cj ~ 0.003; iKennel &: Coronitil Il984al ). where this is not the case. Finally, at these late times 
the spectrum of the PWN discussed in ^ has several features at 7-ray and low frequency {v < 1 
GHz) radio wavelengths which are observed with new facilities such as the Fermi Gamma-ray Space 
Telescope, LOFAR, LWA, MWA, and EVLA. Parameter exploration is required to determine how 
general these conclusions are, but suggest this model will be useful in resolving several outstanding 
problems concerning the evolution of a PWN inside a SNR. 

Additionally, this model is extremely useful in determining if PWNe inside SNRs can de- 
posit a sufficient number of high energy electrons and positrons into the surrounding ISM to ex- 
plain the rising positron f raction of cosmic-ray leptons detected by PAMELA between ~ 1.5 — 100 
GeV (lAdriani et al.ll2009l l and t he excess of cosrnic ray electrons and positrons detected by ATIC 
(|Chang et al.1 bood i and HESS (|Aharonian et ahlbood l between ~ 300 - 800 GeV. These results 
require a ne arby source of high-en ergy electrons and positrons, and one possibility is PWNe inside 



SNRs (e.g. iMalyshev et al.ll2009l ). If this is correct, the average PWN must deposit ~ 10 ergs 



of energetic electrons and positrons into the surrounding ISM, and these particl es must have an 
energy spectrum flatter than E~'^ which extends up to an energy of ~ 1 TeV ( Malyshev et al 
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20091 ). For the set of parameters modeled in ^ these conditions are met for only a short period 
of time during the evolution of this PWN. Using this model, it is possible to determine what sets 
of neutron star, pulsar wind, supernova, and ISM parameters are required for the PWN to satisfy 
these criteria for a longer period of time, and evaluate different models for particle escape from the 
PWN and their effect on the PWN's evolution - particularly if these particles escape gradually or 
suddenly from the PWN. 

One major limitation of this model is that is in inherently one-dimensional since it assumes a 
constant pressure and magnetic field strength inside the PWN. Therefore, it is insensitive to any 
pressure and magnetic field variations inside the PWN. Recent Chandra observations of PWNe 
have revealed the exi stence of internal structure (e.g., the torus and jets in the Crab Nebula ; 



Weisskopf et al.l 12000) but of spectral changes as well (jWeisskopf et al.l 12000; iMori et al.l l2004l ) 
- indicative of a non-uniform pressure and/or magnetic field inside the PWN. We are also in- 
sensitive to possibly significant effects asymmetries in the PWN resulting from either the space 
velocity of the neutron star or inhomogeneities in the progenitor supernova and/or surrounding 
ISM which can have on its evolution, particularly after it collides with the SNR reverse shock (e.g. 



van der Swaluw et al.l 12004 ). Additionally, while we estimate the growth rate of Rayleigh- Taylor 
instabilities in the shell of swept-up material surrounding the PWN, we can not determine the role 
they might play in enabling particles to escape from the PWN - critical in determining if such 
PWNe are responsible for the PAMELA, ATIC, and HESS results discussed above. 

To summarize, we present a one-dimensional model for the evolution of a PWN inside a 
SNR. This model self-consistently evolves the magnetic, dynamical, and radiative properties of 
the PWN throughout its evolution, and therefore represents a significant improvement over other 
currently existing models. The model described here provides a framework for investigating the 
effect of more complicated descriptions of pulsar winds, e.g. the presence of ions in the pulsar 
wind, on the dynamical and radiative evolution of the PWN, as well as determining if there exists 
any unique observational signatures of such processes. Additionally, it is well-suited for using the 
observed properties of a PWN to constrain the properties of the central neutron star, its pulsar 
wind, progenitor supernova, and surrounding ISM. Its applicability to a large number of known, 
well-studied PWNe make it is an extremely powerful for studying these intriguing systems. 
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Zaldariagga for useful discussions. JDG is supported by an NSF Astronomy and Astrophysics 
Postdoctoral Fellowship under award AST-0702957. 
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Table 1: Model Input Parameters 



Parameter 


Units 


Description 






Supernova Explosion 




]^q51 gj,gg 


Initial Kinetic Energy of Supernova Ejecta 




Mo 


Mass oi Supernova Ejecta 






Interstellar Medium 






Numl)or dcuisily of surrounding ISM 






Neutron Star 


p 




Neutron Star Braking Index 


Tsd 


years 


Characteristic Spin-down Timescale of the Neutron Star 


-^0,40 


]^q40 gj,gg g-i 


Initial Spin-down Luminosity of the Neutron Star 


'^psr 


km 


Space velocity of the Neutron Star 






Pulsar Wind 


Ve 




Fraction of neutron's star spin-down luminosity injected as electrons 


m 




Fraction of neutron's star spin-down luminosity injected as ions 


VB 




Fraction of neutron's star spin-down luminosity injected as magnetic energy 


E 


varies 


Minimum Energy of Electrons injected into the PWN 


-E'e,max 


varies 


Maximum Energy of Electrons injected into the PWN 


-^i,min 


varies 


Minimum Energy of Ions injected into the PWN 


^i,mEix 


varies 


Maximum Energy of Ions injected into the PWN 


7e 




Electron injection index 


li 




Ion injection index 
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Table 2: Trial Model Input Parameters 



Parameter 


Value 




1 


Mej 


8 




0.1 


P 


3.00 


Tsd 


500 




1 




120 




0.999 


m 







0.001 


E ■ 


511 keV 


-^e,max 


500 TeV 


-^i,min 








7e 


1.6 


7i 





Table 3: Selected Observed Properties of the Crab Nebula 



Observed Property 


Value 


Reference 


Model Prediction 




1.5-2 pc 


Green f2006l 


1.7 pc 


rts 


0.07-0.14 pc 


WeisskoDf et al. f2000) 


0.24 pc 




~ 1270 km s"^ 


Temim et al. (2006) 


2000 km s"^ 


Radio Luminosity 


1.8 X 10^^ ergs s~^ 


Prail & Scharrinehausen fl997) 


2.2 X 10^^ ergs s"^ 


Q^radio 


-0.3 


Green r2006^ 


-0.3 


-^X,0.5-10keV 


1.3 X 10^^ ergs s~^ 


Mori et al. f2004) 


~ 3 X 10^^ ergs s~^ 


To.S-lOkeV 


1.9-3.0 


Mori et al. f2004l 


1.8 


-^50GeV-50TeV 


~ 10^"^ - 10^^ ergs s"^ 


Aharonian et al. (2004) 


1.6 X 10^^ ergs s"^ 
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Fig. 1. — The ratio of the inverse Compton cross-section of an electron scattering off the CMB 
relative to the Thomson cross-section. 

Fig. 2. — Radius of the SNR (i?snr, blue line), reverse shock {Rrs, green line), PWN (i?pwin red 
line), the position of the neutron star (rpsj., orange line), and the radius of the termination shock 
inside the PWN (rts, purple line) for a neutron star, supernova, and PWN with the properties 
listed in Table [2j The termination shock only exists when the pulsar is inside the PWN, and is 
always centered on the pulsar. The vertical dotted lines indicate the division between the four 
evolutionary phases of the PWN discussed in ^ 

Fig. 3. — The pressure inside the PWN (Ppwm I'ed) and the pressure inside the SNR just outside 
the PWN (-Psnr(-Rpwn)) blue). At times t < tcoi, we assume that -Psnr(-Rpwn) = 0- The vertical 
dotted lines indicate the division between the four evolutionary phases of the PWN discussed in 

m 



Fig. 4. — The velocity of the shell of swept-up material surrounding the PWN (fpwn) red) and the 
velocity of the supernova ejecta just beyond this mass shell [vej{Rpwn), blue]. The vertical dotted 
lines indicate the division between the four evolutionary phases of the PWN discussed in ^ 

Fig. 5. — The mass of material swept-up by the PWN (Msw,pwn) as it expands inside the SNR. 
The vertical dotted lines indicate the division between the four evolutionary phases of the PWN 
discussed in ^ 

Fig. 6. — The spin-down luminosity of the pulsar {Lpsr = E, blue) and the adiabatic (Lad, green), 
synchrotron (Lgynch; red), and inverse Compton (Lie, orange) luminosity of the PWN. Solid lines 
indicate energy lost by the PWN, while dashed lines indicate energy gained by the PWN. The ver- 
tical dotted lines indicate the division between the four evolutionary phases of the PWN discussed 
in 33 

Fig. 7. — The total energy (L^pwm red), particle energy (L/p^n,p5 green) and magnetic energy (L/pwn,B? 
blue) of the PWN. The vertical dotted lines indicate the division between the four evolutionary 
phases of the PWN discussed in ^ 
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Fig. 8. — Magnetic field strength of tlie PWN, i?pwn- The vertical dotted lines indicate the division 
between the four evolutionary phases of the PWN discussed in ^ 



Fig. 9. — The magnetization parameter of the PWN a for t < tco\ (top) and t > tco\ (bottom). 
In the bottom plot, the vertical lines represent the division between the second, third, and fourth 
evolution phase of the PWN discussed in ^ 



Fig. 10. — The energy spectrum of electrons and positrons inside the PWN during the first (upper 
left), second (upper right), third (bottom left), and fourth (bottom right) evolutionary phase discussed 
in ^ In all four plots, the different color lines correspond to the energy spectrum at different ages, 
the vertical dotted lines indicate Ee,min (right) and Ee,ma.x (Wt), and the dashed line at the bottom 
of the plot indicates the shape and extent of the injected spectrum. In the upper-right plot, the 
black line corresponds to the energy spectrum at the time the pulsar leaves its PWN for the first 
time. In the bottom left plot, the dotted orange line corresponds to the energy spectrum at the time 
when the pulsar re-enters the PWN. In the bottom right plot, the dotted yellow line corresponds to 
the energy spectrum just before the pulsar re-enters the PWN. 



Fig. 11. — The spectrum of photons radiated by the PWN during the first (upper left), second 
(upper right), third (bottom left), and fourth (bottom right) evolutionary phase discussed in ^ In 
all four plots, the different color lines correspond to the photon spectrum at different ages, and the 
cross-hatched regions indicate, from left to right, the radio, mid-infrared, optical, soft X-ray, hard 
X-ray, 7-ray, and TeV 7-ray regimes of the electromagnetic spectrum. In the upper-right plot, the 
black line corresponds to the photon spectrum at the time the pulsar leaves its PWN for the first 
time. In the bottom left plot, the dotted orange line corresponds to the photon spectrum at the time 
when the pulsar re-enters the PWN. In the bottom right plot, the dotted yellow line corresponds to 
the photon spectrum just before the pulsar re-enters the PWN. 
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Fig. 12. — The photon luminosity of the PWN in the radio {u = 10^ — 10^^ Hz, Lradioj red), 
mid-infrared (A = 3.6 — 160/xm, orange), near-infrared/optical (A = 2.35/im-354.3 nm, yellow), soft 
X-ray (/iz^ = 0.5 — 10 keV, dark green), hard X-ray {hi' =15 keV-10 MeV, blue), 7-ray {hv = 10 
MeV-100 GeV, dark blue), and TeV 7-ray {hu = 50 GeV - 50 TeV, purple) for t < tcoi (top) and 
t > tco\ ( bottom) . The definition of the wavebands were chosen to reflect the frequency / wavelength 
/ energy coverage of current observing facilities, and are given in units usually associated with that 
portion of the electromagnetic spectrum. In the bottom plot, the vertical dotted lines demarcate 
the second, third, and fourth evolutionary phase discussed in ^ 



Fig. 13. — The spectral index a (photon index F) of the radio (red), 0.5-10 keV (green), and 
2-10 keV (blue) emission from the PWN. The dotted vertical lines indicate the evolutionary phases 
discussed in ^ 



Fig. 14. — The luminosity density of the PWN between the radio and soft X-ray bands during the 
first {upper left), second {upper right), third {bottom left), and fourth {bottom right) evolutionary 
phase discussed in ^ The filled and open stars indicate t'b and L^,^ derived by joining power law 
fits to the radio and 0.5-10 keV (closed stars) or 2-10 keV (open stars) spectrum of the PWN. 



Fig. 15. — The change in spectral index Aa between the radio and 0.5-10 keV band (green) and 
the radio and 2-10 keV band (blue). The dotted vertical lines indicate the evolutionary phases 
discussed in ^ 



Fig. 16. — The total energy inside the PWN (-Epwru red) and the minimum energy inside the PWN 
calculated using the radio and 0.5-10 keV spectrum of the PWN (-E'pwn,min(o.5-iokeV; blue) and the 
radio and 2-10 keV spectrum of the PWN (-Epwn,min(o.5-iokcV) green). The dotted vertical lines 
indicate the evolutionary phases discussed in ^ 



Fig. 17. — Magnetic field strength of the PWN, -Bpwn (red), as well as the magnetic field strength 
inferred using the radio and 0.5-10 keV spectrum of the PWN (i?pwn,min(o.5-iokeV)5 blue) and the 
radio and 2-10 keV spectrum of the PWN (-Bpwn,min(2-iokcV) > green). The dotted vertical lines 
indicate the evolutionary phases discussed in ^ 
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Fig. 18.— The minimiiiii cingula,r scale unstcible to Rayleigli-Taylor instctbilities (^rt,crit ~ 
Art,crit/-Rpwn; red) and the angular scale maximally unstable to Rayleigh- Taylor instabilities 
(^rt.max = Art,max/-Rpwn; bluc) for t < icol (top) and t > tcol (bottom). In the bottom plot, the 
vertical dotted lines demarcate the second, third, and fourth evolutionary phase discussed in ^ 



Fig. 19. — The growth rate {uj^t, Equationl3ip of Rayleigh- Taylor instabilities with different angular 
scales for t < icoi (top) and t > icoi (bottom). In the bottom plot, the vertical dotted lines demarcate 
the second, third, and fourth evolutionary phase discussed in ^ 



Fig. 20. — The penetration depth of Rayleigh- Taylor bubbles into the PWN {h^t, orange) and radius 
of the PWN (i?pwni red). The dotted vertical lines indicate the evolutionary phases discussed in 

m 
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A. Evolution of a Non-Radiative Supernova Remnant 



In this model, we assume that the progenitor core-collapse supernova injected material with 
mass Mej and initial kinetic energy Egn into a constant density pisM interstellar medium (ISM) 
environment. We assume the supernova ejecta initially comprise a constant density inner core 
surro unded by a p oc outer envelope - the standard assumption for core-collapse supernovae 
(e.g. iBlondin et al.l l200lh - wher e the transition velocity vt between these two components is 



(Equation 3 in 



Blondin et al 



2001 



Vt 



18Mei 



1/2 



Under this assump tion, the density of material inside the SNR evolves as (Equations 1 
Blondin et ;il]l200ll ): 



Pei{r, t) 



10 p 



10 



Vtt 



r < Vtt 
r > Vtt 



(Al) 



2 in 



(A2) 



Initially, the evolution of the SNR is self-similar (e.g. IChevalieij Il982l ) . and therefore described in 
terms of characteristic length {Rch), time (ich)) and mass { Mrh) scales. For a SNR e xpanding into 
a constant density ISM, these scales are (Equations 1-3 in lTruelove h McKedll999l ): 



(A3) 
(A4) 
(A5) 



The expanding supernova ejecta drives a shock wave into the surrounding ISM (called the "forward 
shock") which marks the outer boundary of the SNR (i?snr)- Initially the dynamics of the SNR is 
dominated by expanding supernova ejecta because the ejecta mass is much greater than the mass of 
the ISM material swept u p and shocked by the SN R (Msw,snr = t^-^snrPiSM)- During this period, 
i?snr(i) is (Equation 75 in 



Truelove fc McKeelll999l ^: 



t V"^ 

iisnr(t) = 1.12i?ch(— j 



(A6) 



where 1.12 is specific for a /j oc r ^ outer ejecta envelo pe, though it varies little for different 



values of the power-law exponent ([Truelove &: McKedll999l). Therefore, the expansion velocity v^ 
(vsnrit) = ^%^) of the SNR is fEauation 76 in lTruelove McKeelll999l l: 



0.75 



tch 



t 

tch 



-1/3 



(A7) 



^sw,snr ~ -^ej • 



As the SNR grows, -Mswsnr increases and eventually will reach the point where Afg, 
At this point, the swept-up ISM material will begin to dominate the dynamics of the SNR, and 
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the SNR is said to enter the Sedov- Taylor phase of its evokition (jSedovl Il959l : iTavloij llQSOl ) . For 
a SNR with the ejecta p rofile given in Equation \A2\ this transition occurs at ti me ^st — 0-52^rh 
dXruelove McKeelll999l l. During the Sedov- Taylor phase, i?snr is (Equation 56 in lTruelove &: McKed 
I999I ): 



r5/2 



2.026- 



1/2 



PismJ 



{t - tsT) 



2/5 



(A8) 



where -Rsnr,ST = -Rsnr(isT)- The Sedov- Taylor phase ends when the gas recently s hocked by the 



expanding SNR cools radiatively, which occurs approximately at (Equation 3 in iBlondin et al 
19981 ): 



^rad 



2.9 



4/17 



]^q51 ergs^ 

after the supernova, where rup is the mass of a proton. 



-9/17 



X 10 years 



(A9) 



The pressure P and density p profile of a SNR changes significantly over time. Initially, the 
expansion velocity of the SNR is significantly larger than the sound speed of the surrounding ISM. 
As a result, the swept-up material is shocked by the surrounding ejecta. Assuming that energy 
losses due to cosmic ray acceleration is negligible, the pressure of the recently shocked ISM material 

[Psnr(-Rsnr)] is: 



^ Pism^snr 



it? 



(AlO) 



assuming an adiabatic index 7 = 5/3 (non-relativistic gas). The pressure of the shocked ISM is 
significantly higher than that of the ejecta inside the SNR, which due to adiabatic expansion has 
cooled significantly since the explosion. As a result, the layer of shocked ISM expands inside the 
SNR, driving a shock wave, referred to as the reverse shock, into the supernova ej ecta. When the 
reverse shock is in the o uter envelope of the supernova ejecta, its radius R^s is (jChevalierill982l : 
Truelove fc McKeelll999l ): 



Rrsit) 



1 



1.19 



Rsnv (^) ) 



and therefore, its velocity relative to the surrounding ISM Vj^s is: 



.it) 



1 



1.19 



rit). 



(All) 



(A12) 



The properties of ejecta recently shocked by the reverse shock depends on the velocity of the 
reverse shock relative to the unshocked ejecta, v^y The standard assumption is the unshocked 
ejecta is expanding ballistic ly \vpj(r,t) = r/t]. There fore, the relative velocity of the reverse shock, 
Vrs = Vej{Rrs,t) - Vrsit), is (jTruelove k McKeelll999l ): 



2.38 



rit). 



(A13) 
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Eventually, the reverse shock will enter the constant density core at the center of the SNR. For a 
SNR with an ini tial ejecta density profile d efined in Equation (jA2|), this occurs a time icoro — 0.25tch 
(Equation 79 in lTruelove &: McKed 11993 ) where tch is defined above in Equation (IA4D. After this 
time, the radius of the reverse shock evolves as (Equation 83 in lTruelove &: McKeelll999l ): 



i?rs(i) 



1.49-0.16- 



.t-t. 



tch 



0.46 In 



^ch 
tch 



t. 



During this stage, v^s is equal to (Equation 84 in lTruelove &: McKeelll999l ): 

, (t ^core 



0.50 + 0.16- 



^ch 



^ch 



(A14) 



(A15) 



Since the velocity of the reverse shock is much greater than the sound speed in the unshocked 
ejecta, the reverse shock - like the forward shock - is a strong shock. Ther efore, the pressure of 
ejecta recently shocked by the reverse shock Prsi is (jTruelove &: McKeelll999l ): 

PcjiRrs,t)Vj.sit)'^ 



Prs{t) 



"-fsnr {t^ . 



(A16) 



PlSMWsnr(i)^ 

When the reverse shock is in the outer envelope of the supernova ejecta, the density, velocity, 
and pressure profile of the SNR between the reverse sh o ck an d forward shock (i?rs < r < iJsnr) 
is calculated using the equations derived by IChevalieii (jl982l ). These equations are no longer 
valid once the reverse shock enters the constant density ejecta core. At this point, we model the 
density, velocity, and press ure structure of the SNR using the solution for a Sedov- Taylor SNR 
presented in Appendix A of iBandieral (|1984| ). This is valid until the SNR goes radiative {t = trad)- 
Unfortunately, the internal structure of a radiative SNR is poorly understood, and therefore we 
do not attempt to model the evolution of a PWN in this environment. Ho wever, previous work 
suggests that the radius of a PWN .Rpwn inside a radiative SNR evolves as (IBlondin et al.ll200ll ; 
van der Swahiw fc WJI2OO1I ): 



pwn ^ ^0.075 



(A17) 



B. Initial Properties of a PWN inside a SNR 

In order to model the evolution of a PWN inside a SNR, it is necessary to estimate the initial 
conditions. To estimate the initial energy of the PWN, we assume that adiabatic losses dominate 
and the PWN is expanding with a constant velocity. In this case: 

^ = i^ofl + -V^-^. (Bl) 
at V '^sd/ t 

Defining y = {p + l)/(p — 1), e = -£'pwn/(^o'^sd)) and x = t/Tgd, one derives that: 

^ = -i + (i + x)-y. (B2) 

ax X 
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Using the boundary conditions that e(x = 0) = (initiahy, the PWN has zero energy), one derives 
that, for y = 2 (p = 3): 

ln(l + 2;) 1 



x + 1 



and, for y / 2 (p / 3), 



(1 + x) 



1-y 



+ 



x{l-y)i2-y) x(l-y)(2-2/)- 



(B3) 



(B4) 



We also assume that initiahy the fraction of the PWN's energy in magnetic fields is r/^, the fraction 
of the PWN's energy in electrons and positrons is rje, and the fraction of the PWN's energy in ions 
is r]i. Additionally, we assume that initial spectrum of particles inside the PWN has the same shape 
as the injection spectrum. 

If adiabatic losses dominate, the equation of motion of the PWN for Psnr(-Rpwn) ^ is 
(IQstriker Gunnlll97ll : Ichevalier fc Franssonlll992l ): 



dv 



pwn 



swjpwn 



dt 



^^^pwn[-^pwn PcjiRpwii) ^ (^pwn ^cj(-^pwn)) ]; 



where: 



dE, 



pwn 



dt 



E 



-fpwn^TT^p^jj tipwn • 



(B5) 



(B6) 



At early times {t <C r), E k, Eq, and it is possible to solve Equations IB5I fc |B6I analvticallv (e.g. 
Chevalier &: Franssonlll992l ). For the init ial supernova ejecta density given in Equa tion I A2 1 Rp^„ 



can b e expressed as (e.g. Equation 2.6 in lChevalier &: Franssonlll992l . Equation 5 in lBlondin et al 

2nnih: 



i?p„„(t) = 1.44 
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